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We consider 0(a s ) corrections to the squared masses of the pseudo-Goldstone excitations about 
the ground state of dense quark matter. We show that these contributions tend to destabilize the 
vacuum, leading to a surprisingly complex phase structure for quark matter as a function of quark 
mass, even for small a s . In particular we find two new phases of CFL quark matter possibly relevant 
for the real world, for which #qcd = tt/2. 
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The dependence of the QCD vacuum on the masses of 
the light quarks is most efficiently analyzed by means of 
a chiral Lagrangian, which allows one to systematically 
compute corrections to the chirally symmetric vacuum in 
a power expansion in the quark masses. The utility of the 
chiral Lagrangian approach is that it properly accounts 
for the lightest QCD excitations, the pseudo Goldstone 
bosons (PGBs), whose masses vanish in the chiral limit 
and which therefore have the most important role in de- 
termining the vacuum structure for small quark mass. 
In addition to allowing one to compute perturbative cor- 
rections to the vacuum structure, the chiral Lagrangian 
also allows one to investigate phase transitions at critical 
values of the quark masses. An example of such a phase 
transition was first discovered by Dashen pj, who found 
a line of phase transitions in the (m«,mrf) plane corre- 
sponding to spontaneously broken CP symmetry, where 
m Uj d are the up and down quark masses respectively. 

In hadronic matter at nonzero baryon density, more 
complicated phase transitions have been investigated us- 
ing chiral Lagrangians, including both pion |2| and kaon 
condensation 0, Q. More recently, chiral Lagrangians 
have been used to analyze the vacuum structure of dense 
quark matter near the chiral limit. It has been convinc- 
ingly argued that for degenerate quarks such a system 
should be color superconducting y|, and that for three 
massless flavors the color-flavor-locked (CFL) ground 
state is preferred 0, 0, 0- Chiral perturbation theory 
can be used to describe the ground state and excitations 
of CFL matter away from the chiral limit 0, 0, [TT| . 
It is already known that the phase structure is much 
more complicated than in the vacuum case. In partic- 
ular, an analysis to leading order in both quark masses 
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and the QCD coupling constant a s reveals a phase tran- 
sition to a kaon condensed phase along a line in the m s -m 
plane 000, where m s is the strange quark mass, 
and for simplicity we are considering the isospin limit, 
m u = m c i = m. In this Letter we show that the leading 
order calculation does not capture the full complexity of 
the quark matter ground state, and that there exist new 
phases revealed at sub-leading order in a s at small quark 
mass in the m s -m plane which could be relevant for the 
real world. 

The reason that 0(a s ) effects are important is that 
the leading order result exhibits an accidental degener- 
acy which is only lifted at next-to-leading order. Once 
0(a s ) effects are included the meson masses receive neg- 
ative corrections to their squared masses which can be 
larger than the leading order contributions, even for small 
a s . In particular, a meson of mass M with M 2 cx m 2 
at leading order can become destabilized by a pertur- 
bative correction of the form SM 2 cx ~a s m 2 sl given that 
m 2 3> m 2 . Characteristic of these new phases, which can 
coexist with kaon condensation, is that they break an ex- 
act Z2 symmetry (parity) and that det U ^ 1, where U is 
the unitary matrix characterizing the quark condensate. 
Such phases appear to have a non-vanishing strong CP 
violating angle, Qqcd = tt/2. 

In the CFL phase, attractive interactions near the 
(shared) Fermi surface cause quarks of all three flavors 
to undergo BCS-like pairing. The resulting condensate 
takes the form 
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where A3.g are the pairing gaps in the antisymmetric 
color triplet and symmetric color sextet channels respec- 
tively. A perturbative calculation yields 



Ai = a 



(In 2) 

1627T 



-A 



(2) 



The gap A3 is larger than A6 because one-gluon exchange 
is attractive in the triplet channel but repulsive in the 



sextet channel; the radiative generation of nonzero Ag 
does not affect the symmetries of the CFL phase. The 
condensate spontaneously breaks the U(1)b baryon num- 
ber symmetry, as well as the SU(3) C x SU(3)l x SU(3)r 
symmetry down to the diagonal subgroup, SU(3) c +l+r 
where SU(3) C refers to the color gauge group. As a result 
the gluons become massive, and there appears a nonet of 
Goldstone bosons analogous to the pseudoscalar meson 
nonet in the vacuum (in addition to the superfluid mode 
from the breaking of U(1)b)- Unlike in the vacuum, in- 
stantons are suppressed at high density, and the entire 
nonet is expected to be light. 

In the CFL phase low energy excitations may there- 
fore be parametrized by E = e 2%lx l^ , V = e 4 "' , and 

B = e iP/fB witn n = n a T a^ the T a bemg 577(3) gen _ 

erators. The n a and rj' form the pseudo-scalar nonet, 
while j3 (which we shall ignore for the remainder of this 
Letter) is the Goldstone boson of broken baryon num- 
ber. Under the original symmetry of SU (3)l x SU (3)r x 
U(1)b x U{1)a the E field transforms as (3,3W ), V 
transforms as (1, l)(o,4)j where we have assigned axial 
charges of +1 and —1 to the right- and left-handed 
quarks respectively. A gauge invariant order parameter 
for the system is (schematically) (qlQlQrQr) whose ori- 
entation in SU(3) L x SU(3)r x U(1)a/SU(3) v is given 
by the U(3) matrix U = VE. The quark mass matrix 
M = diag(m, m, m s ) explicitly breaks the chiral symme- 
tries of QCD, and enters the effective theory as a spurion 
transforming as (3,3W_ 2 ) © (3, 3)( 0i2 )- 

The chiral Lagrangian introduced in [j| and elaborated 
in [Tot [TTL fl3| is an expansion in quark mass, derivatives, 
and strong coupling a s . Generically it takes the form 
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where A = A 3 . To date the theory has been discussed 
expanding the dimensionless function C to include all 
terms of order (<9 2 /A 2 ), (m 2 /^A) 2 , and (M 2 /u 2 ); the 
Lagrangian to this order is 

C = ^Tr (DoEA^ - w 2 TrVE • Ve) 
f 2 ' 

+ -7^ (d VdoV^ - v 2 VV ■ VV*) 

+ a 3 V f [(TrME) 2 - TrMEME] + h.c. (4) 

The covariant derivatives in Eq. J3J contain the Bedaque- 
Schafer term ^| that acts as a dynamically generated 
chemical potential for flavor symmetries, 



L> E = <9 E + i [M 2 /2u, E] 



(5) 



By making use of the effective theory introduced by Hong 
[l5| , all of the parameters in Eq. @ have been computed 
to zeroth order in a s [Tfl HH 13 • The results we will use 
in this paper are 
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FIG. 1: Feynman diagrams corresponding to the g M terms 
in the effective potential. The open and solid squares de- 
note insertions of the gap, and the dot is an effective vertex 
which describes chirality changing quark-quark scattering am- 
plitudes. 



The a 3 term in Eq. which gives rise to the me- 
son masses, can be rewritten as 2a 3 U TrME + h.c. where 
M = diag(mm s , mm s , m 2 ). It is therefore apparent that 
there exists a meson with square mass proportional to 
m 2 A 2 //i 2 , corresponding to fluctuation in the {3, 3} com- 
ponent of VE, which is exceedingly light. A major point 
of this paper is that it does not make sense to keep a term 
proportional to m 2 while dropping ones proportional to 
a s m 2 , given that (m/m s ) 2 ~ 1 x 10~ 3 . There are various 
small parameters that control the chiral expansion in the 
CFL phase. Defining (m 2 //^A 3 ) = 5 and m 2 //! 2 = e, 
we take m s m/fj, 2 ~ 5 and a s m 2 /[i 2 ~ e with S 2 <C e <C 6, 
expanding the Lagrangian to 0(S) and 0(e). 

Fig. 1 shows the leading order graphs that contribute 
to the 0(M 2 ) terms in the chiral effective lagrangian. 
Fermion lines in the graph on the left correspond to the 
anomalous quasi-particle propagators in the CFL phase 
and depend on the gaps A 3 6- The four- fermion vertex 
corresponds to a chirality changing four-fermion interac- 
tion of order 0(g 2 M 2 ) The Feynman diagram on 
the left in Fig. 1 gives rise to the a 3 term in in the low 
energy Lagrangian Eq. , as well as the new term 



8Cq = -a 6 V f [(TrME) 2 + TrMEME] + h.c. , (7) 
a 6 = A( A e) 2 , (8) 



which is proportional to a s , given the relation Eq. 
[l7j . Perturbative corrections to the Feynman diagram 
in Fig. 1 that are not summed by the gap equation con- 
tribute at 0(a s mm s A 2 ). We neglect these terms in this 
work. We also neglect the contribution from the graph on 
the right of Fig. 1 which is of order 0(?n 2 A 4 //i 2 ), which 
is exponentially small as compared to the terms we have 
kept. 

The contributions to the vacuum energy proportional 
to a 3 and ae are opposite in sign, reflecting the fact that 
one gluon exchange is attractive in the color triplet chan- 
nel, and repulsive in the sextet channel. Because of the 
minus sign in SCq one finds that the 0(a s m 2 ) contribu- 
tions to the square of the meson masses are negative, and 
can alter the vacuum alignment. 

To analyze the vacuum alignment, we consider free en- 
ergy density of the homogeneous medium arising from 
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Neglecting the new contributions, one finds that f2 
is minimized by the vacuum configuration 



V = 1 
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cos# isinf 
1 i sin 9 cos 9 , 



(10) 



where the angle 9 is satisfies 



cos# = min [1, cq] 



co 



16ma3 



(11) 



where we have dropped 0(m/m s ) corrections to Co- This 
is the kaon condensed phase discussed in 0, 0, El- 
Note that if we neglect the a§ contribution, then since 
^33 = there exists a very light excitation if about 
T,k, corresponding to 



VT, 



J K 
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or 



, g-iy/3 

V = e lv/3 , E = I e-^/ 3 
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The ip excitation has a mass M 2 oc {a^m 2 / p, 2 ) in the ab- 
sence of kaon condensation, which is very small. There- 
fore it is in this direction that we expect an instability to 
arise when the clq operator is included, which contributes 
5M 2 ~ — (asm 2 / fj, 2 ). With the ansatz Eq. (|T?S|> . we can 
expand the free energy to first order in a s and m/m s , 
with fi = f2o + Hi + • • - j where 



Oo = - 



! 2 m\ sin 2 9 



— 4mm s 03 cos 
2 „ 2 , 



8p 2 

fii = —4 (03m 2 cos 9 — a & m 2 cos 2 9) cos y> . (14) 



We are neglecting contributions to the free energy of or- 
der mm s ae, m 2 m 2 and smaller. 

Minimizing fig with respect to 9 determines the K° 
condensate. The subleading contribution to the free en- 
ergy, Ox, has a negligible effect on the kaon condensate, 
which is still given by equ. Qll|). but is the leading contri- 
bution to the potential for the angle tp. Minimizing with 
respect to tp yields 



cos ip = Sign [a^m 2 — a^m 2 cos#] 



(15) 



m 




m Q 



FIG. 2: The phase diagram as a function of light quark mass 
m and strange quark mass m s . The phases marked CFL, K° , 
77, and r/ + K° are respectively the CFL phase without meson 
condensation, with kaon condensation, with r\ condensation, 
and with both r\ and K° condensation. Phase transitions are 
represented by a solid line if first order, a dashed line if second 
order. The location of the tetracritical point is {m* s ,m*}, 
given in Eq. Ilbl : including subleading effects causes it to 
split into two tricritical points. 



For brevity we will refer to the phase with cos<p = — 1 
as an "77" condensate, although the mode that condenses 
is actually a linear combination of 77 and if . (Note that 
in this phase V = e i7r / 3 ; since V has U(1)a charge equal 
to four, this phase is equivalent to a strong CP angle 
8qcd — Tt/2.) Thus we see that we have four different 
phases, depending on the quark masses, corresponding 
to whether or not the kaon and/or the 77 condense. 

The symmetry structure of the four phases is as fol- 
lows. The ordinary CFL phase (neglecting isospin break- 
ing) possesses an H = SU(2) x U(1)y x P symmetry, 
where P is parity, under which E — > E^ and V — > V'. 
Neutral kaon condensation by itself breaks H down to 
J7(l)em x P' 1 where P' is a combination of ordinary par- 
ity and a discrete hypercharge rotation exp(iirY), and 
?7(l)em is the electromagnetic group. The 77 condensate 
by itself breaks the discrete parity symmetry P, but no 
continuous symmetries |l8j . The phase with both 77 and 
K° condensation breaks H to U(l) em . 

The phase diagram is easily constructed by means of 
Eqs. (|11I15|) . The phase boundary for kaon condensa- 
tion is second order and is described by an equation of 
the form m oc m 3 ; it is unaffected by the existence or 
absence of 77 condensation at the order in a s and m/m s 
that we are working. The phase boundary for 77 conden- 
sation is first order, as one might expect for the breaking 
of a discrete symmetry. However, to the order we work 
the ip angle does not encounter a free energy barrier at 
the phase transition, and so it is possible that higher or- 
der corrections could render the phase transition second 
order. It is described by a line m oc m s in the region 
without kaon condensation, and the curve m oc l/m s in 
the kaon condensed region. The phase diagram shown in 
Fig. [3 At this order there is a tetra-critical point, which 
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is separated into two tricritical points when higher order 
corrections are included. The coordinates of this critical 
point are (m*,m*) given by 



2(ln2) 
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If the baryon chemical potential is very large we can use 
perturbation theory in order to estimate the location of 
the tetracritical point. If we assume that there are only 
three flavors of quark, and that a s (m T ) — 0.353 as in the 
real world, then at (i = 10 10 GeV we have a s /ir — 0.009, 
and the critical values 



H = 10 10 GeV 



= 320 keV 77i* = 89MeV . (17) 



In the regime /_( ~ 0.5 GeV, relevant to the physics of neu- 
tron stars, the coupling is large and perturbation theory 
is not reliable. If we assume that a s /ir ~ 1, A ~ 100 
MeV, and f„ ~ 90 MeV we find 



500 MeV : 



4.6 MeV 



120 MeV 



(18) 



It is interesting to note that the values of the critical 
masses are remarkably close to the physical values, par- 
ticularly in the strong coupling example. As a conse- 
quence, we expect that the region of the phase diagram 
explored in this work and displayed in Fig. 2 may well be 
relevant to the structure of CFL matter in neutron stars. 

We emphasize that if CFL quark matter is the true 
ground state of quark matter at large baryon density 
then our study provides a rigorous analysis of the phase 
diagram for small quark masses in the regime m < m s . 
As such, it provides the correct starting point for an ex- 
ploration of the phase structure for larger values of the 
quark masses, and an important guide for any attempt 
to numerically simulate QCD in the regime of very large 
baryon density and low temperature. 
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